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Objectives 


Damage  detection  in  engine  bladed  disks  is  often  performed  through  ultrasound  and  eddy 
current  techniques  that  are  reliable,  but  expensive  and  lack  in-situ  monitoring  capability. 
Alternatively,  vibration-based  damage  detection  methods  are  relatively  inexpensive,  have  real¬ 
time  in-situ  potential,  but  are  generally  inaccurate  due  to  low  sensitivity.  The  goal  of  this 
research  is  to  advance  the  state-of-the-art  of  vibration-based  damage  detection  of  bladed  disks  by 
utilizing  the  unique  vibration  localization  characteristics  of  such  periodic  structures  to  enhance 
damage  detection  sensitivity  and  robustness  through  piezoelectric  circuitry  networking. 


Summary  of  Efforts 

In  this  research,  we  have  successfully  developed  an  innovative  methodology  to  enhance 
sensitivity  and  robustness  of  vibration-based  damage  detection  for  periodic  structures  such  as 
engine  bladed  disks.  A  new  piezoelectric  network  concept  to  temporarily  increase  the  vibration 
localization  in  periodic  structures  has  been  established.  Analytical  and  experimental  studies 
have  been  carried  out  to  obtain  insight  towards  the  underlying  principles  of  the  proposed 
approach.  Fundamental  understandings  of  vibration  energy  propagation/distribution  in  periodic 
systems  with  and  without  piezoelectric  circuitry  are  reached  via  wave  and  vibration  analyses. 
Design  guidelines  of  circuitry  elements  that  can  optimize  the  system  performance  have  been 
established.  Finally,  multivariate  statistical  analysis  tools  are  synthesized  to  evaluate  and 
demonstrate  the  improvement  in  sensitivity  and  robustness  of  damage  detection  owing  to  the 
localization  enhancement  of  the  piezoelectric  networking. 


Descriptions  of  Accomplishments  and  New  Findings 

The  timely  detection  of  damage  in  spatially  periodic  structures  such  as  engine  bladed  disks  is 
an  extremely  challenging  task.  While  ultrasound  and  eddy  current  technologies  have  good 
reliability,  they  need  significant  amount  of  human  involvement,  have  narrow  field  coverage,  and 
are  position  sensitive.  When  applied  to  engine  systems,  they  generally  require  the  complete 
disassembly  of  engines,  which  leads  to  very  long  maintenance  time  and  high  cost.  On  the  other 
hand,  the  advent  of  the  blade-tip-timing  (BTT)  sensing  technology  is  an  important  progress  in 
turbomachinery  bladed  disk  sensing.  This  technology  uses  a  number  of  probes  installed  in  the 
engine  casing  to  sense  the  time  instants  at  which  the  blades  pass  the  probes.  When  properly 
analyzed,  these  blade  passing  times  can  lead  to  the  vibratory  amplitude  information  at  the  blade 
tips.  There  has  been  recent  exploration  in  using  BTT  sensing  for  vibration  response-based, 
online  bladed  disk  damage  detection. 

While  vibration-based  approaches  such  as  those  using  frequency-shift  information  to  infer 
damage  occurrence  have  been  successfully  utilized  in  many  applications,  periodic  structures  such 
as  bladed  disks  have  clustered  natural  frequencies,  which  makes  it  difficult  to  ascertain  the 
change  of  these  frequencies  caused  by  damage.  Another  problem  that  is  related  to  the  detection 
performance  is  the  uncertainty  and  variation  of  the  structures.  In  particular,  bladed  disks,  albeit 
manufactured  with  ultra-high  precision,  all  possess  certain  level  of  manufacturing  and 
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assembling  error.  The  blade-to-blade  random  difference  is  often  referred  to  as  mistuning.  If  one 
uses  the  frequency-shift  information  alone,  potentially  fatal  damage  could  be  masked  by  the 
inherent  mistuning  in  a  bladed  disk. 

When  the  coupling  between  the  substructures  is  weak  in  a  periodic  structure,  it  is  well  known 
that  mistuning  may  induce  the  so-called  vibration  localization  in  the  structure,  where  a  single  or 
several  substructures  experience  excessively  large  vibration  response.  For  ideally  periodic 
structures  without  mistuning,  the  substructural  coupling  will  split  the  identical  substructural 
natural  frequencies  into  groups  of  frequencies  and,  in  any  given  vibration  mode,  all  the 
substructures  will  have  the  same  response  amplitude  but  with  certain  phase  difference.  These 
frequency  groups  are  referred  to  as  passbands.  The  frequency  ranges  outside  the  passbands  are 
the  stopbands.  Within  a  frequency  stopband,  wave  propagation  throughout  the  structure  exhibits 
spatial  attenuation.  The  width  of  a  passband  is  associated  with  the  level  of  coupling.  In  a 
weakly  coupled  situation,  the  passband  width  is  small,  i.e.,  the  modal  density  is  high.  When  the 
system  has  mistuning,  the  passband  behavior  is  affected.  Narrower  passband  is  generally  more 
susceptible  of  losing  the  passband  characteristics  due  to  mistuning.  Essentially,  in  such  situation, 
the  reflected  portion  of  the  spatially  propagating  wave  due  to  mistuning  discontinuity  between 
the  substructures  is  more  significant,  and  the  successive  reflections  of  the  wave  propagation  lead 
to  vibration  localization  for  frequencies  inside  the  original  passband. 

Obviously,  uncontrolled  vibration  localization  results  in  local  concentration  of  vibration 
energy,  and  thus  is  harmful  to  the  system  durability.  Nevertheless,  some  features  of  localized 
vibration  in  a  periodic  structure  can  potentially  benefit  damage  detection.  Weak  coupling  and 
high  modal  density  make  the  system  highly  sensitive  to  mistuning  including  the  damage  effect. 
In  such  case,  even  though  the  natural  frequency-shift  caused  by  damage  is  still  small  and  may 
not  be  noticeable  when  inherent  mistuning  is  present,  the  vibration  response  patterns  will  be  very 
different  upon  damage  occurrence.  In  other  words,  for  a  system  with  strong  vibration 
localization,  even  a  small  damage  superimposed  to  the  inherent  normal  mistuning  in  the  structure 
may  cause  the  vibration  pattern  to  be  significantly  different  from  that  of  the  healthy  structure. 

The  overall  goal  of  this  research  is  to  fundamentally  advance  the  statc-of-thc-art  of  damage 
detection  of  periodic  structures,  by  developing  a  highly  sensitive  and  robust  approach  exploiting 
the  unique  vibration  localization  potential  of  such  structures.  Specifically,  an  innovative 
piezoelectric  circuitry  networking  methodology  that  can  temporarily  induce  or  intensify 
structural  vibration  localization  to  amplify  the  damage  effect  on  the  system  vibratory  signature 
during  the  inspection  stage  is  created.  The  network  can  be  easily  switched  off  (open  circuit) 
under  normal  operating  conditions.  The  overall  concept  is  to  form  an  electro-mechanical  wave 
channel  to  enhance  localization  in  damaged  structures  through  energy  redistribution  and  wave 
reflection.  Under  the  temporarily  intensified  localization,  the  system  response  becomes 
extremely  sensitive  to  damage  occurrence  (which  will  further  enable  robust  detection),  since 
even  a  small  damage  can  now  cause  drastic  difference  in  vibratory  response  patterns.  This 
amplified  feature  change  is  then  analyzed  using  a  suite  of  robust  decision  making  algorithms  that 
can  deal  with  measurement  noise/uncertainty.  With  this  new  innovation,  a  vibration-response- 
based  health  monitoring  system  can  become  effective  for  periodic  structures  such  as  bladed  disks, 
entertaining  the  merits  of  being  simple  and  cost  efficient,  and  can  be  applied  to  online 
monitoring  when  combined  with  effective  sensing  techniques. 
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During  the  course  of  this  research,  systematic  research  activities  have  been  carried  out, 
including  experimental  analysis,  analytical  characterization,  and  numerical  investigation.  The 
following  steps  and  milestones  are  accomplished: 

■  Evaluation  of  the  underlying  principle  and  development  of  tools  for  integrated  system 
synthesis; 

■  Creation  of  guidelines  for  network  parameter  design  to  maximize  vibration  localization  to 
amplify  the  damage  effect; 

■  Examination  and  demonstration  of  damage  detection  performance  with  variation  and  noise. 
The  research  results  and  new  findings  are  highlighted  in  the  following  sections. 


Part  I.  Underlying  principle  and  feasibility  study 

We  consider  a  system  shown  in  Figure  1.  For  illustration  and  without  loss  of  generality,  the 
periodic  structure  is  assumed  to  consist  of  N  identical  cantilever  beams/blades  coupled  with  N 
springs.  Here  the  springs  are  employed  to  emulate  the  mechanical  coupling  effect  in  generic 
periodic  structures,  such  as  the  stiffness  coupling  due  to  disk  in  an  engine  bladed  disk.  To 
induce  or  enhance  temporarily  the  vibration  localization  effect  during  inspection,  identical 
piezoelectric  inductive 
circuits  are  integrated  onto 
all  substructures/blades, 
where  in  each  circuit  an 
inductor  is  connected  in 
series  with  the  piezoelectric 
transducer  embedded  to  the 
substructure.  These  circuits 
are  coupled  through 
identical  capacitance 

elements  to  form  a  global  Figure  1  (a)  Schematic  of  periodic  structure;  (b)  Periodic 

network  (Figure  lb).  structure  integrated  with  piezoelectric  network. 


Our  interest  is  in  the  global  dynamic  response  of  the  periodic  structure.  While  each 
substructure  has  an  infinite  number  of  degrees-of-freedom,  we  employ  the  assumed  mode 
method  to  describe  the  substructure  motion.  For  simple  illustration,  here  we  only  use  the  first 
substructure  mode  for  discretization.  Let  ^  be  the  first  local  beam  mode  without  the 
piezoelectric  circuit.  The  transversal  displacement  of  the  7-th  blade  is  approximated  as 
Wj{xJ)  =  (^{x)qj(t),  where  qj{t)  is  the  generalized  mechanical  displacement  (y  =  The 

equations  of  motion  of  electro-mechanically  integrated  system  can  be  obtained  as 


+  kqj  +  l^pqQj  +  K  (^j  -  9y-i )  +  K  -  9j+\ )  =  fj 

Qj-i  - Qj^\ ) + :r ^  ^pp^j  ^ ® 


(la) 

(lb) 


For  the  original  mechanical  structure  without  the  piezoelectric  network,  the  equation  of  motion  is 
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mijj+gqj+kqj  +K{qj-qj.i)  +  k^{qj  =  fj 


(2) 


In  the  above  equations,  m ,  g  and  k  are  the  equivalent  mass,  damping  and  stiffness  of  the 
mechanical  substructure,  is  the  electrical  charge  flow  in  the  7-th  circuit,  kpp  is  the  inverse 
capacitance  of  the  piezoelectric  transducer,  kp^  indicates  the  electro-mechanical  coupling,  L  is 
the  circuitry  inductance,  and  fj  is  the  external  disturbance  on  the  y-th  blade.  We  define 
k^  =  1/C^  as  the  inverse  of  the  coupling  capacitance.  Under  harmonic  excitation  at  frequency 
CO,  the  force  and  displacements  can  be  expressed  as  fj  =  fje''"',  qj  =qje'‘^'  and  Qj  =Qje'‘“’ . 
Here  /  =  . 

The  above  equations  describe  the  perfectly  periodic  system  without  mistuning.  In  reality, 
however,  there  are  always  imperfections  especially  in  the  mechanical  structure.  As  the  common 
treatment  in  localization  study,  we  assume  that  the  structural  mistuning  only  exists  in  the 
substructure  stiffness.  The  stiffness  of  the  /-th  blade  with  mistuning  can  then  be  expressed  as 

kj  =  k  +  5k^j,  where  Sk^j  is  the  random  stiffness  mistuning  with  zero-mean.  The  damage  is 
assumed  to  occur  on  the  /-th  blade  only,  causing  stiffness  loss  Sk^i,  and  hence  the  stiffness  of  the 
/-th  blade  with  damage  is  kj  =  k  +  Sk^j  +  Sk^/.  For  non-dimensionalization,  wc  define  v  =  coj 

where  (o„  =  y[k/m  .  The  equations  of  motion  of  the  system  with  and  without  the  piezoelectric 
circuitry  become,  respectively. 

-n^x  +  2g(ni)Xj  +  (1  +  ) X  +  (x  -  x  (x  -  x  )  +  <^Syj  =  —p-f=  (3a) 

co„^m 

-n\yj  +  rI  {2yj  -  yj_\  -  (2^y  -  -  ^^>1 )]  +  ^Sxj  +  S'^y^  =  0  (3b) 

and 

-Cl^Xj  +  2g{Qi)Xj  +  (1  +  hsj  )xj  +  R]  {x^  -  )  +  R^  {x^  -  )  =  0  (4) 

where 

=V^pp/'^’  ^  ,  S  =  0)J(0„,  ^  -  kp^l  ,Jk^ ,  Xj  =  yfmqj,  yj  =  yfiQj 

Q.  =  col(o„,  R^  =  yli^,  R^  =  yjka/kpp  ,  As^  =Ak„j/k 

Here  co„  is  the  original  mechanical  resonant  frequency  of  the  blade,  co^  is  the  natural  frequency 
of  the  electrical  circuit  g  is  the  structural  damping  ratio,  S  is  essentially  the  inductance  tuning 
ratio.  ^  is  the  generalized  electro-mechanical  coupling  coefficient  which  reflects  the  energy 


5 


transfer  capability  of  the  piezoelectric  transducer,  and  yj  are  the  generalized  mechanical  and 
electrical  displacements,  respectively,  Q  is  the  non-dimensional ized  excitation  frequency;  is 
the  non-dimensionalized  coupling  between  the  blades,  and  is  the  coupling  capacitance  ratio. 
ASj  is  the  mistuning  ratio  which  is  a  zero-mean  random  number  with  stander  deviation  cr . 

For  the  /-th  blade  with  damage,  the  corresponding  equations  are 

-Cl^x,  +  2g(Qi)x,  +{\  +  As,+Ak, )x,  +  (x,  -  x,_, )  + R^{x,-x,^^)  +  4Sy,  =  /'  (5) 

co-y/m 

-Q^Xi  +  2g{Q.i)Xi  +(1  + A.S7  +  +  Rc(^i  ®  (^) 


where  Akj  is  the  stiffness  loss  ratio,  -1  <  Akj  <  0  . 

The  original  mechanical  substructures  (beams/blades)  all  have  one  degrec-of-freedom.  With 
the  substructural  coupling,  in  an  ideal  situation  without  mistuning,  the  identical  blade  resonant 
frequencies  split  into  a  group  of  natural  frequencies  of  the  periodic  structure.  This  frequency 
group  forms  a  frequency  passband.  From  the  wave  propagation  perspective,  passband  is  the 
frequency  range  within  which  the  wave  propagates  throughout  the  substructures  without  spatial 
attenuation.  When  mistuning  is  present,  the  passband  behavior  is  affected.  When  coupling  is 
weak  and  the  original  passband  is  narrow,  mistuning  will  more  easily  eliminate  the  passband 
characteristics  and  lead  to  more  severe  vibration  localization.  Generally,  vibration  localization  is 
due  to  the  high  sensitivity  of  vibratory  response  with  respect  to  mistuning  including  damage. 
Therefore,  our  underlying  principle  is  that,  if  we  can  temporarily  increase  the  level  of  vibration 
localization  during  inspection,  we  can  highlight  damage  occurrence. 


We  have  carried  out 
experimental  investigation  to 

demonstrate  this  underlying 

principle.  The  overall  test  setup 
with  a  periodic  structure 

emulating  bladed-disk  is  shown  in 
Figure  2.  The  structure  tested  is  a 
6-bay  specimen.  On  each 
substructure,  we  placed 

piezoelectric  transducers  to  act  as 
exciters  to  emulate  engine  order 
excitation.  The  blades  are 

obviously  randomly  mistuned.  To 
emulate  irregularity  such  as 

damage,  we  added  a  small  mass  to 
one  of  the  blades.  The  results  are 
shown  in  Figure  3.  We  have 

examined  cases  that  are  weak  and 


Figure  2  Experimental  setup  for  underlying  principle  examination. 
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strong  in  localization  to  start 
with.  For  weakly  localized  case, 
there  are  very  small  differences 
before  and  after  damage  occurs, 
and  the  differences  are  generally 
in  the  noise  level  and  are  not 
detectable.  However,  for  the  case 
that  we  have  strong  localization 
to  start  with,  the  healthy  and 
damaged  structural  responses  are 
very  different,  as  clearly 
demonstrated  in  Figure  3.  This 
exp)eriment  has  clearly  verified 
our  underlying  principle  that 
localization  enhancement  can 
greatly  improve  the  vibration- 
based  damage  detection 
performance. 

Part  II.  Design  of  the  piezoelectric  circuitry  network 

In  Part  I  of  this  research,  we  have  demonstrated  the  underlying  idea  of  vibration  localization 
enhancement  to  benefit  damage  detection.  For  the  original  mechanical  structure,  if  the 
mechanical  coupling  is  not  weak,  the  system  is  less  sensitive  to  the  mistuning  and  damage,  and 
the  vibration  localization  is  less  significant  even  though  local  wave  reflection  still  occurs.  The 
integration  of  piezoelectric  inductive  circuit  onto  the  substructure  adds  an  additional  degree-of- 
freedom  to  each  substructure.  After  these  circuits  are  coupled  through  capacitance  elements 
(Figure  lb),  we  obtain  two  frequency  passbands,  one  corresponding  to  the  original  mechanical 
frequencies,  and  the  other  corresponding  to  the  circuitry  dynamics  with  electro-mechanical 
coupling.  The  widths  and  locations  of  these  two  frequency  passbands  are  determined  by  a 
number  of  parameters  including  the  circuitry  elements.  Our  hypothesis  is  that,  by  properly 
designing  the  piezoelectric  circuitry  network,  we  can  obtain  a  narrow  frequency  passband  that  is 
much  more  sensitive  to  substructural  change  (i.e.,  damage)  compared  to  the  original  mechanical 
frequency  passband.  In  other  words,  for  frequencies  within  this  narrow  passband,  the  vibration 
localization  will  be  intensified. 

Wave  propagation  analysis  and  passband  bandwidth 

The  integrated  system  has  two  wave  propagation  channels,  one  through  mechanical  coupling 
between  the  substructures,  and  the  other  through  electrical  coupling  in  the  piezoelectric  network. 
For  each  coupling  type,  there  will  be  a  right-going  wave  component  and  a  left-going  wave 
component.  For  the  perfectly  periodic  structure,  within  the  frequency  passbands,  the  waves 
propagate  freely  with  no  spatial  attenuation.  The  natural  frequencies  of  the  periodic  structure  arc 
all  inside  the  passbands.  Outside  the  passbands,  the  propagating  wave  amplitude  would  decay 
e.xponentially.  One  effective  way  to  analyze  the  passband  feature  is  the  wave  transfer  approach. 
For  the  integrated  system,  we  define  the  displacement  vector  of  the  y-th  interface  between  the 
adjacent  subsystems  as  Uy  Xj  yy]^,  which  consists  of  both  the  mechanical  and 
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Normalized  blade  amplitude  w/  localization 


Figure  3  Experimental  result  of  maximum  normalized  blade 
amplitude  change  due  to  damage. 


electrical  generalized  coordinates  from  the  y-th  and  the  (/+l)-th  subsystem.  Neglecting  the 
damping,  we  can  re-write  Equation  (3)  in  a  transfer  matrix  format,  =T^u^_, ,  where  is  a 

4x4  transfer  matrix  given  as 


1  +  2RI  + 

^3 


n^R 


2d2  g 


1  +  Ay^  -Q 


2  A 


yj 


1 
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1  ^ 
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2  +  -V 

R^ 

a  J 

^  / 

0 
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-1 

0 

0 

0 


0 


0 

0 


(7) 


Obviously,  the  transfer  matrix  is  a  function  of  excitation  frequency.  For  the  mistuned  system, 
the  transfer  matrix  is  random  due  to  the  random  mistuning  ratio  Asj.  For  the  ideally  periodic 

system  without  mistuning,  the  transfer  matrix  is  identical  for  all  subsystems.  For  the  original 
structure  without  the  piezoelectric  network,  we  have  v^  =S^Vy_, ,  where  \j  =  Xj^  ,  and 

the  transfer  matrix  is  given  as 


\  +  2R]+^SJ-Q^ 


-1 

0 


(8) 


Previous  investigations  on  periodic  structures  have  used  the  Lyapunov  exponent  to  analyze 
the  passband  features  and  vibration  localization.  Mathematically,  the  Lyapunov  exponent  is 

obtained  as 


Non-dimensional  frequency 
Figure  4  Lyapunov  exponent  of  original  mechanical 
structure  without  mistuning,  Rc  =  0.1. 


y(uo)=  lim -^logfu^l  (9) 

Here  Ug  is  the  displacement  vector  of  the  first 
substructure,  and  's  the  displacement  vector 
of  the  A-th  substructure.  The  number  of 
Lyapunov  exponents  is  equal  to  twice  the 
number  of  the  coupling  coordinates.  For  the 
original  mechanical  system  with  mechanical 
coupling  only,  there  are  two  Lyapunov 
e.xponents  which  come  in  pair.  For  the 
integrated  system  analyzed  in  this  research,  there 
are  totally  four  Lyapunov  exponents  which  come 
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in  pairs.  The  positive  Lyapunov  exponents  describe  the  exponential  decay  of  vibration 
amplitudes.  Inside  the  passbands,  the  Lyapunov  exponent  is  zero,  indicating  no  decay.  We  then 
use  the  Wolf  s  algorithm  to  calculate  the  Lyapunov  exponents  of  the  periodic  system. 

We  first  analyze  the  original  mechanical  system  without  mistuning.  The  generalized  electro¬ 
mechanical  coupling  coefficient  is  ^  =  0.2  and  the  inductance  ratio  is  selected  as  ^  =  1  (this 
selection  will  be  explained  in  the  next  section).  Here  we  assume  the  structure  has  strong 
mechanical  coupling,  i.e.,  i?^  =  0.1.  As  illustrated  in  Figure  4,  there  exists  one  frequency 

passband.  We  then  analyze  the  system  integrated  with  the  piezoelectric  network.  In  order  to 
introduce  weak  electrical  coupling  into  the  integrated  system,  the  capacitance  value  is  selected  as 
=0.0913,  which  will  be  further  discussed  in  the  analysis  that  follows.  As  shown  in  Figure  5, 

the  integrated  system  with  piezoelectric  network  has  two  frequency  passbands.  owing  to  the 
introduction  of  the  piezoelectric  network. 


Non-dimensional  frequency  Non-dimensional  frequency 

Figure  5  (a)  Lyapunov  exponent  of  system  integrated  with  piezoelectric  network,  without 
mistuning;  (b)  Zoomed-in  I**  passband;  (c)  Zoomed-in  2^  passband. 
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Figure  6  Eigenvalue  change  with  different  Ra  (6=1)  ■  -  :  smallest  one  in  the  mode  group; 

—  —  •  :  largest  one  in  the  mode  group,  (a)  1"*  mode  group;  (b)  2”**  mode  group 


As  mentioned  in  Part  1,  the  passband  width  generally 
reflects  the  sensitivity  of  the  integrated  system  with 
respect  to  mistuning  or  damage  occurrence.  The  width 
of  the  passband  is  closely  related  to  the  electrical 
elements  in  the  piezoelectric  circuitry.  For  tuned  system 
without  mistuning,  the  edges  of  a  specific  passband  are 
determined  by  the  smallest  (marked  by  solid  line)  and 
the  largest  eigenvalues  (marked  by  dashed  line)  in  the 
corresponding  modal  group,  plotted  in  Figure  6.  The 
eigenvalues  change  with  varying  capacitance  ratio  . 
We  can  see  that:  (1)  for  each  passband,  its  bandwidth 
reaches  its  minimal  value  at  certain  ;  and  (2)  the 
minimal  passband  width  occurs  in  the  2"^  passband 
under  R^  =0.0913,  which  is  referred  to  as  the  optimal 
value  (Figure  7). 


x10‘^ 


respect  to  Ra. 


Inductance  tuning  effect  and  forced  response  magnitude 

The  actual  damage  detection  will  be  carried  out  using  forced  responses  that  can  be  measured 
by  the  BTT  measurement.  The  excitation  form  studied  is  the  engine  order  excitation  that  is 

commonly  used  in  bladed  disk  analysis,  i.e.,  fj  =  .  Here  /q  is  the  force  magnitude, 

(f>j  =  {2n{E  -  \){j  -  X)]/ N  is  the  phase  of  the  force  applied  to  the  y-th  blade,  and  E  is  the  engine 

order  number  {E  =  \,  --,N).  For  a  perfectly  periodic  system,  under  a  certain  engine  order,  only 
the  corresponding  vibration  mode  is  excited.  For  a  mistuned  system  or  a  system  with  damage, 
however,  all  modes  are  excited  even  under  one  specific  engine  order,  and  excessive  response 
amplitude  may  appear  as  a  result  of  vibration  localization.  The  frequency  response  of  all  blades 
will  then  be  used  in  damage  detection. 
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We  now  proceed  to  analyzing  the  inductance  tuning  effect.  Recall  that  the  non-dimensional 
inductance  is  defined  as  S  =  cOglco^  ,  the  ratio  of  the  local  circuit  natural  frequency  to  the  un¬ 
coupled  blade  natural  frequency.  Indeed,  the  inductance  tuning  directly  decides  the  local  circuit 
natural  frequency,  and  therefore  dictates  the  dynamic  coupling  between  the  mechanical  structure 
and  the  piezoelectric  network.  The  forced-response  plot  of  a  representative  blade  (the  5*  blade) 
when  we  select  5  =  0.5  is  shown  in  Figure  8a.  In  this  case,  the  first  resonance  is  mainly  related 
to  the  piezoelectric  circuitry  dynamics  and  is  away  from  the  mechanically  dominant  resonance. 
Here  we  can  see  that  this  additional  resonance  due  to  circuitry  is  20  dB  lower  than  the 
mechanically  dominant  resonance.  When  we  select  ^  =  1 .5 ,  the  second  resonance  is  mainly 
related  to  the  piezoelectric  circuitry  dynamics  and  is  away  from  the  mechanically  dominant 
resonance.  Again,  the  additional  resonance  due  to  circuitry  is  20  dB  lower  than  the  mechanically 
dominant  resonance  (Figure  8b).  When  we  select  ^  =  1 ,  as  the  local  circuit  has  the  same  natural 
frequency  as  the  un-coupled  blade  natural  frequency,  the  two  resonant  peaks  have  similar 
response  magnitudes,  as  shown  in  Figure  8c.  Therefore,  we  conclude  that  the  strong  dynamic 
coupling  at  <5  =  1  ensures  that  the  additional  resonant  responses  due  to  circuitry  can  be 
realistically  measured,  which  is  used  as  the  optimal  inductance  tuning  throughout  this  research. 


Excitation  Frequency  (HZ) 


Figure  8  Resonant  peak  amplitude  comparison  under  different  inductance  tunings,  (a)  5=0.5; 

(b)  6=1.5;  (c)  6=1. 


II 


Localization  enhancement  with  proper  capacitance  limins 

The  preceding  analyses  on  passband  width  and  response  measurements  are  carried  out  under 
the  condition  that  the  system  is  ideally  periodic  without  mistuning  or  damage.  We  now 
investigate  the  effect  of  mistuning  including  damage  under  the  optimal  inductance  and 
capacitance  tunings  (i.e.,  S  =  1  and  =0.0913).  Figure  9b  plots  the  Lyapunov  exponents  of 

the  nearly  periodic  structure  (with  mistuning  standard  deviation  cr  =  0.001)  with  different 
coupling  capacitance  ratio  /?„  within  the  original  frequency  passband.  For  comparison  purpose, 
the  Lyapunov  exponents  of  the  ideally  periodic  structure  are  plotted  in  Figure  9a.  The  non-zero 
Lyapunov  exponents  indicate  that  the  corresponding  wave  propagation  will  have  amplitude 
decay,  thereby  causing  vibration  localization.  Different  values  affect  the  Lyapunov 


e.xponents  differently.  One  can  see  that  R^  =0.0913  affects  the  Lyapunov  exponent  the  most. 

leading  to  its  largest  increase,  i.e.,  the  largest  spatial  amplitude  decay  in  wave  propagation.  This 
clearly  demonstrate  that  the  vibration  localization  can  be  maximized  when  we  select 
R^  =  0.09 1 3  that  yields  the  minimal  passband  width. 


Figure  9  Lyapunov  exponents  at  the  original  2"^  passband  frequencies,  (a)  Tuned 
system  without  mistuning;  (b)  Mistuned  system  (a  =  0.001) 

—  •  -  :  Ra  =  0.085;  -  :  Ra  =  0.0913  (optimal  value); - -  :  Ra  =  0.098. 


We  further  analyze  the  vibration  modes  of  the 
systems  to  gain  insight  to  vibration  localization 
enhancement.  We  first  consider  the  original  mechanical 
structure  without  the  piezoelectric  network.  For  the 
structure  with  mistuning  (with  mistuning  standard 
deviation  <7  =  0.001),  the  original  mechanical  coupling 
between  the  substructures  7?^  =0.1  is  quite  large.  We 
also  assume  that  damage  occurs  on  the  20*  blade  with 
stiffness  loss  Ak2Q  =  0.2%  .  We  randomly  pick  the  5* 

mode  for  demonstration  and  plot  it  in  Figure  10.  Here 
the  dash-dot  line  is  for  the  perfectly  periodic  structure; 
the  dashed  line  is  for  the  mistuned  structure  without 
damage,  and  the  solid  line  is  for  the  mistuned  structure 
with  damage.  For  the  perfectly  periodic  structure,  the 
mode  is  spatially  harmonic,  extended  throughout  the 


Blade  number 

Figure  10  S’"  Vibration  mode  of  mechanical 
structure  without  piezoelectric  network. 

—  •  —  :  Perfectiy  periodic  structure; 

- :  Mistuned  periodic  structure; 

-  ;  Mistuned  structure  with  damage. 


12 


structure.  When  mistuning  and  damage  occur,  the  mode  is  no  longer  spatially  harmonic. 
However,  since  the  mechanical  coupling  is  not  weak  in  this  case  (./?^  =0.1),  the  mode  of  the 

mistuned  structure  does  not  exhibit  strong  localization.  Furthermore,  no  obvious  pattern 
difference  is  observed  after  the  damage  occurrence. 

We  then  analyze  the  system  integrated  with  the  piezoelectric  network  (  J  =  1  and 
=0.0913).  As  mentioned,  the  integrated  system  has  two  modal  groups.  Figure  1  la  shows 

the  5*’’  mode  for  the  tuned  system,  mistuned  system  and  damaged  system.  This  S'**  mode  belongs 
to  the  first  modal  group  and  the  first  passband.  It  is  obvious  that  no  severe  localization  can  be 
observed  under  mistuning,  and  no  apparent  mode-shape  pattern  change  exhibits  when  damage 
occurs.  The  reason  is  that  under  the  tuning  5  =  1  and  =0.0913,  the  first  passband  width  is 

not  minimized  (Figure  6a).  Figure  lib  shows  the  45'*'  mode,  belonging  to  the  second  modal 
group  for  the  tuned,  mistuned  and  damaged  system.  It  is  worth  emphasizing  that  this  modal 
group/passband  reaches  the  minimal  width  under  the  tuning  ^  =  1  and  R^  =0.0913  (Figure  6b 

and  Figure  7).  Clearly,  severe  vibration  localization  can  be  observed  under  mistuning,  and  the 
mode-shape  pattern  has  significant  change  when  damage  occurs.  These  results  are  consistent 
with  the  preceding  Lyapunov  exponent  analysis  and  wave  propagation  analysis,  and  demonstrate 
the  effect  of  piezoelectric  circuitry  on  localization  enhancement.  With  these  analyses,  the 
optimal  design  of  the  piezoelectric  circuitry  has  been  accomplished  and  validated. 
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Figure  11  Vibration  mode  of  structure  integrated  with  piezoeiectric  network,  (a)  S®"  mode  from 
the  1*'  modai  group;  (b)  45®'  mode  from  the  2"®  modal  group. 

—  •  —  :  Perfectly  periodic  structure; - :  Mistuned  periodic  structure;  - - -  :  Mistuned 

structure  with  damage. 

Part  III.  Damage  detection  performance  investigation  and  robust  decision  making 

In  this  part  of  research,  we  examine  actual  damage  detection  performance  by  using  blade  tip 
responses  obtained  during  engine  spin-up  or  spin-down  processes,  through  simulated  numerical 
investigations.  For  the  damage  detection  purpose,  we  observe  the  system  response  within  a 
certain  frequency  range.  Since  generally  the  damage  effect  is  most  significant  around  the 
resonances,  we  analyze  the  response  magnitudes  under  frequency  sweep  around  the  resonant 


13 


frequencies.  While  mistuning  is  present  in  all  blades,  damage  only  occurs  on  one  blade.  In  order 
to  characterize  the  forced  response  difference,  we  define  a  performance  index  which  is  the 
normalized  ratio  of  the  difference  between  the  forced  responses  before  and  after  damage 
occurrence  within  certain  frequency  range  (from  to  CO2)  around  the  resonances.  For  the y-th 


Figure  12  Damage  index  definition. 

:  frequency  response  curve  of 
healthy  baseline; 

- - :  new  measurement  to  be 

examined. 


blade,  the  mathematical  form  of  the  index  called 
“Normalized  Relative  Difference"'  (NRD)  is  expressed  as 

^  C  -  "i )  ( '  0) 

As  shown  in  Figure  1 2,  the  dashed  line  is  the  frequency 
response  curve  before  damage,  and  the  solid  line  is  the  one 
after  damage;  is  the  shaded  area  between  curves 

and  Xj^  .  Clearly,  the  higher  this  index  value,  the  more 

significant  the  change  in  vibration  pattern  (due  to  damage) 
for  the y-th  blade  is. 


Performance  improvenient  due  to  localization  enhancement  via  circuitry  integration 

In  this  illustration,  we  let  =  40 ,  i.e.,  the  number  of  blades  is  40.  The  non-dimensionalized 
mechanical  coupling  ratio  is  selected  as  =0.l .  The  inherent  mistuning  ratio  ASj  has  standard 

deviation  cr  =  0.00l  .  The  damping  ratio  is  selected  as  f  =  0.001  .  Compared  with  the 
mechanical  coupling  ratio,  the  mistuning  level  is  low,  which  implies  that  the  original  mechanical 
structure  does  not  exhibit  significant  vibration  localization. 
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Figure  13  Normalized  damage  Index 
comparison  when  k2o=0.1% 

Original  structure  w/o  piezo  network 
H  :  Integrated  structure  with  piezo  network 


Figure  14  Normalized  damage  index 
comparison  when  k2o=0.3% 

Original  structure  w/o  piezo  network 
H  :  Integrated  structure  with  piezo  network 


A  small  damage  on  the  20^  blade  is  assumed  to  causes  stiffness  loss  AA:2o  =  0.1%.  We 

begin  our  analysis  with  the  mechanical  structure  without  the  piezoelectric  circuitry,  and  excite 
the  structure  with  the  third  engine  order  excitation  to  obtain  the  forced  response  of  each  blade. 
With  the  response  curves  before  and  after  damage  for  the  20^  blade,  we  can  calculate  the  values 
ofNRDj  (j  =  1,*»*,40)  for  all  blades.  The  NRD  values,  plotted  in  Figure  13  for  the  mechanical 
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system  without  the  piezoelectric  network,  are  all  below  1 .38.  The  results  corresponding  to 
stiffness  loss  M20  =0.3%  are  plotted  in  Figure  14.  For  this  more  severe  damage,  the  NRD 
values  for  the  mechanical  system  without  the  piezoelectric  circuitry  network  are  all  below  4.12. 


We  then  study  the  response  of  the  bladed-disk  integrated  with  the  piezoelectric  network 
under  the  aforementioned  optimal  circuitry  tuning  (^  =  1  and  =0.0913).  We  again  analyze 

the  forced  response  of  each  blade  under  the  third  engine  order  excitation  with  frequency  sweep. 
The  integration  of  the  piezoelectric  circuitry  introduces  an  additional  degree-of-freedom  to  each 
substructure.  The  integrated  subsystem  now  has  two  resonances.  We  focus  our  attention  on  the 
second  resonance.  The  reason  for  choosing  this  resonance  to  analyze  is,  as  explained  in  Part  II, 
this  second  resonance  is  within  the  passband  with  minimized  width  that  has  the  highest  modal 
sensitivity.  We  calculate  NRDj  (y  =  I, •••,40)  for  the  aforementioned  two  damage  cases  and  plot 
the  results  in  Figure  13  and  Figure  14.  In  both  cases,  with  the  integration  of  the  piezoelectric 
network,  the  overall  NRD  values  are  significantly  increased.  The  maximum  is  3.28  for 
Mjo  =0.1%  and  9.06  for  M20  =  0.3% .  The  more  severe  the  damage,  the  larger  the  difference  in 

the  NRD  value  is.  The  comparisons  in  Figure  13  and  Figure  14  have  demonstrated  that  the 
integration  of  piezoelectric  circuitry  network  has  more  than  doubled  the  maximum  damage  index. 


Statistical  analysis  and  robust  decision  makine  under  variation  and  noise 

The  preceding  analysis  demonstrates  the  performance  improvement  under  one  specific  set  of 
mistuning.  Since  the  mistuning  is  random  in  nature,  we  carry  out  Monte  Carlo  simulation  on 
blade  mistuning  to  investigate  the  performance  from  a  statistical  perspective.  Figure  1 5  plots  the 
maximum  damage  index  under  Monte  Carlo  simulation,  when  M20  =  0.1% .  The  horizontal  axis 


is  the  maximum  damage  index,  and  the  vertical  axis  is  the  result  occurrence  frequency  under 
Monte  Carlo  sampling.  Without  the  piezoelectric  network,  the  mean  result  of  maximum  damage 
index  is  1.16,  and  the  standard  deviation  is  0.09.  Now  with  the  integration  of  the  piezoelectric 
network,  the  mean  is  3.27,  and  the  standard  deviation  is  0.06.  Clearly,  the  mean  of  the  maximum 
damage  index  mean  is  increased  nearly  three  times,  and  the  detection  results  have  smaller 
standard  deviation.  These  both  indicate  significant  performance  improvement. 


Maximum  damage  index 

Figure  15  Monte  Cario  simuiation  resuit  of  maximum  damage  index  under 

random  mistuning 
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In  reality,  the  timely  detection  of  damage  is  always  complicated  by  operating  uncertainty 
such  as  response  measurement  noise.  Therefore,  we  have  taken  the  measurement  noise  into  the 
analysis,  and  synthesized  a  series  of  multivariate  statistical  analysis  tools  to  achieve  robust 
decision  making.  We  employ  the  Principal  Component  Analysis  (PCA)  technique  for  data 
compression,  feature  extraction  and  de-noising.  We  first  obtain  a  group  of  response 
measurement  from  the  healthy  bladed-disk  system  to  set  up  the  baseline.  We  then  implement 
PCA,  which  is  the  eigen-analysis  of  variance  on  baseline,  to  separate  feature  space  from  noise 
effect.  Once  a  new  measurement  is  obtained,  we  can  de-noise  it  by  eigen-truncation  to  truncate 
the  noise  content,  and  calculate  the  damage  indices.  Finally,  confidence-level  based  fault 
detection  is  facilitated  by  employing  the  Hotelling  statistical  analysis,  which  is  two-step 
process.  Firstly,  we  evaluate  the  statistics  of  the  baseline  healthy  information  and  set  up  a 
control  limit  with  confidence-level,  e.g.,  99%.  In  step  two,  we  calculate  the  statistic  of 
damage  indices  with  new  measurement  and  compare  with  the  control  limit.  Damage  occurrence 
will  be  declared  with  confidence-level  if  the  statistics  of  certain  damage  indices  exceeds  the 
control  limit  specified. 

Figure  16  is  the  result  comparison  under  damage  =  0.1%.  The  measurement  has  1% 

noise  with  zero  mean.  Figure  16a  is  for  the  result  without  the  piezoelectric  network,  and  the  T‘ 
statistics  of  all  damage  indices  are  below  the  control  limit,  which  means  we  cannot  declare 
damage  occurrence  under  such  noise  level.  Figure  16b  is  for  the  result  after  we  integrate  the 
piezoelectric  network,  and  the  T^  statistics  clearly  exceeds  the  control  limit,  which  indicates 
damage  occurrence  with  99%  confidence  level.  Therefore,  under  the  same  noise-level,  damage 
can  be  successfully  detected  with  the  integration  of  the  piezoelectric  network.  This  demonstrates 
that  the  new  scheme  can  greatly  increase  the  robustness  of  damage  detection. 
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Figure  16  Robust  decision  making  using  T2  statistical  analysis  (a)  without  piezoelectric  network;  (b)  with 
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